Analysis of the resistance in p-SiGe over a wide temperature range 
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The temperature dependence of a system exhibiting a 'metal-insulator transition in two dimensions 
at zero magnetic field' (MIT) is studied up to 90K. Using a classical scattering model we are able to 
simulate the non-monotonic temperature dependence of the resistivity in the metallic high density 
regime. We show that the temperature dependence arises from a complex interplay of metallic and 
insulating contributions contained in the calculation of the scattering rate 1/td(E,T), each domi- 
nating in a limited temperature range. 
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To date the microscopic origin of the metallic phase of 
what has become known as the 'metal-insulator transi- 
tion in two dimensions at zero magnetic field' (MIT) is 
still under discussion Q. We have recently analyzed the 
temperature dependence of the metallic phase in p-SiGe 
H in terms of temperature dependent screening j|]|] 
and quantum corrections to the conductivity ||. At low 
temperatures the temperature dependence of the conduc- 
tivity could be quantitatively described by three contri- 
butions, namely a(T) = <Jry{T) + 5<Jwl{T) + 8<7t(T), 
where ct_d(T) is the Drude conductivity, 8<Jwl{T) and 
5(Jr(T) are the weak localisation (WL) and the inter- 
action contributions, respectively. It was shown that 
the metallic behavior in the SiGe-systems stems from 
screening effects contributing to <jd{T). In the present 
study we focus on the classical Drude part of the con- 
ductivity and extend the experimental range of temper- 
atures up to 90 K. We aim at reproducing the non- 
monotonic temperature dependence of the experimental 
curves with a classical scattering description of the re- 
sistivity y. In this way we are able to make a quanti- 
tative comparison with the experimental data in a wide 
range of temperatures (1.7K < T < 20K) and densities 
(1.5 x 10 n cm~ 2 < p < 5.1 x 10 n cm- 2 ). Below 4K tem- 
perature dependent screening is the dominant contribu- 
tion to the temperature dependence of the resistivity. At 
higher temperatures the smearing of the Fermi function 
and phonon scattering produce a non-monotonic temper- 
ature dependence. The good agreement of theory with 
experiment indicates that temperature dependent screen- 
ing cannot only account for the metallic behavior close 
to the metal-insulator transition but is equally important 
in adjacent temperature and density regimes. The model 
cannot be applied in the insulating regime (fc^D < 1, 
where k-p is Fermi wave vector and Id Drude mean free 
path). 

In our MBE-grown samples the 2DHG resides in a 20nm 
Sio.85Geo.15 quantum well sandwiched between two un- 
doped Si layers. Remote Boron doping was introduced at 
a distance of 15nm above the quantum well and a Ti/Al 



Schottky gate allowed tuning the hole density. For the 
transport measurements standard Hall-bars were used 
(L=300/im, W=300^m). The experiments were carried 
out in a 4 He-system in a temperature range of 1.7K - 
90K, using standard four terminal AC-technique. The 
hole density in the 2DHG was adjustable from complete 
depletion to a density of 5.1xlO n cTO -2 , yielding a maxi- 
mum mobility of 7800cm 2 /Vs at the highest density. The 
ratio of the Drude scattering time td and Shubnikov- 
de Haas relaxation time r q is close to one, indicating 
that large angle scattering dominates the transport in 
such structures ||Q . The insulating contributions of the 
weak localization effect and electron-electron interaction 
become significant at temperatures < IK. For a detailed 
study of these effects see Rcf . 0] . 
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FIG. 1. Temperature sweeps for different densities 

p - 5. 1-2. Ox lO^cm -2 ; the long range sweep (90K-1.7K, thin 
line) is done at B=0.5T, the short range sweep (10K-1.7K, 
thick line) at B=0T. The WL-contribution acts localizing and 
is suppressed at B=0.5T. • represent the temperature T c 
above which the system becomes non-degenerate, ■ marks 
the Dingle temperature Td- 



Fig. p] shows measurements of the longitudinal resis- 
tance p xx as a function of temperature between 1.7K 
and 90K for different densities p. A typical curve (e.g. 



p = 2.3xl0 11 cm 2 , see Fig. g) can be divided into four 
different temperature ranges: 

i. T<4K: the temperature dependence is metallic. 

ii. 4K<T<14K: the temperature dependence is non- 
monotonic and a turnover to insulating behavior 
at higher temperatures occurs. 

iii. 14K<T<60K: the dependence changes again to 
metallic behavior. 

iv. 60K<T<90K: in this regime the 2DHG is not the 
only conducting channel in the sample (visible in 
the Hall density, not shown) . This leads to a strong 
resistivity decrease as the temperature is raised 
from 60K to 90K . 
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FIG. 2. 
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A closeup of the temperature dependence of the 



resistivity for the density p = 2.3x10 cm" . Metallic and 
insulating behavior alternates depending on the temperature 
range. 



In order to achieve a detailed understanding of these 
results we apply the standard semi-classical scattering 
model for the resistivity MM which is valid in the limit 
of /cfZd > 1- 

For finite temperatures and zero magnetic field the resis- 
tivity is expressed as 



Pxx{T) 



pe 2 {r{T)Y 



with the effective mass of m* — 0.25mo (obtained from 
the T-dependence of Shubnikov-de Haas oscillations) and 
the scattering time 



(r(T)) 



i r°° 

— / dEEr(E,T) 
Ef Jo 



df(E,p,T) 
dE 



(1) 



where f(E, \x, T) = (e^-^/^+l)- 1 is the Fermi distri- 
bution function (/i temperature dependent chemical po- 
tential, k-Q Boltzmann constant) and Ep is the Fermi 
energy. 
The transport scattering rate is calculated according to 



1 __ m* 
t{E,T) ~ ^ 



dO 



\V{q E )\ 



e 2 (q E ,T) 



(1 -cos 9), 



(2) 



where q = ^4m*E(l — cos8)/H is the momentum trans- 
fer in a scattering event. 

In the calculation we consider ionized impurity scattering 
H , interface roughness scattering g , alloy disorder scat- 
tering ||, ionized background impurity scattering |1C| ] 
and acoustic phonon scattering pi[ . They enter into 
< |V(g£;)| > which represents the sum of all individ- 
ual mechanisms. 

In this formulation the T-dependence of p xx arises due 
to the energy averaging in Eq. m and due to the explicit 
T-dependence of Lindhard's dielectric function 



e(q, T) = l + V(q)U(q, T)F(q) [I - G(q)} , 



(3) 



where V(q) = 2ire 2 /eq is the interaction potential of the 
hole gas in two dimensions. The T-dependence in Eq. 
again is due to energy averaging when calculating the 
polarizability function [12] 
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U(q,T)= / dE U(q, T = 0, E) 
Jo 



df(E,fi,T) 
dE 



(4) 



In Eq. ||F , ((7) is the form factor for the inversion layer m, 
and G(q) = q/g\/q 2 + k'p with a degeneracy factor g is 
the local field correction in the Hubbard-approximation 
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3. Polarizability function II(g, T) for various temper- 



atures and the density p = 2.3 x 10 cm 

Let us look in more detail into the Eqs. and H. We 
can assert two general effects for the temperature depen- 
dence: (1) the temperature dependence of the chemical 
potential n{T) = k B ,T\ii(e E "/ k ^ T - 1): this shifts the 
balance point around which the integral is performed to 
lower energies with increasing temperature (2) the inte- 
gration itself performs an energy averaging around fJ,(T) 
with width of about 4k%T and depends strongly on the 
curvature of the integrand. In both cases the structure 



is the same: the integrand is a product of the symmet- 
ric (around E — fi) function df j dE with an asymmetric 
function, i.e. Et and II, which in total yields asymmetric 
functions but with opposite sign in the curvature 
The curvature of II(g, T) is negative (see Fig. H) mean- 
ing that dll/dT < =» de/dT < => dr/dT < due 
to energy averaging. The temperature dependence of the 
chemical potential leads to a slight damping of the en- 
ergy averaging effect but can not change its overall be- 
havior. In our case it is negligible: the effect becomes 
only important for temperatures higher than T c at which 
the chemical potential /_(-T) deviates about 10% from its 
zero temperature value (i(0) — E-p (indicated by filled 
circles in Fig. |l|). 
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FIG. 4. The influence of different scattering contributions 
V{q) on the integrand of Eq. m Er e (E,T), for different tem- 
peratures. 

Due to reduced screening t(E,T) in Eq. nj decreases 
with increasing temperature, but energy averaging coun- 
teracts this tendency since the curvature of Et(E,T) is 
positive and therefore d (t) jdT > 0. We illustrate this 
behavior of Et(E,T) in Fig. El At the lowest temper- 
atures the temperature dependence of the polarizability 
function leads to a metallic behavior in the resistivity. In 
the special case of large angle scattering, i.e. q = 2fcp 
which is dominant in SiGe samples jl, the tempera- 
ture dependence of LI(<7, T) is strongest. Only at higher 
temperatures the insulating contribution of the energy 
averaging becomes significant and finally overcomes the 
metallic screening effect. 

We now discuss the interplay of these effects in the tem- 
perature ranges (i)-(iii) mentioned in the beginning: 

i. 2DHG is degenerate: the screening function dom- 
inates the temperature behavior. For large angle 
scattering (i.e. scattering for q r_ 2k f) it has been 
shown that |JHJ] 



a D (T)=a D (0) 1-C(p) 
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counteracts temperature dependent screening. At 
higher temperatures energy averaging dominates, 
i.e. the temperature dependence is insulating. 

iii. the system becomes non-degenerate: phonon scat- 
tering is the strongest scattering mechanism; the 
observed temperature dependence is due to phonon 
freeze out. This behavior is in agreement with 
the generally accepted linear temperature depen- 
dence of acoustic deformation potential scattering 
for which a larger temperature coefficient of the re- 
sistivity is expected for lower carrier densities g. 

Deviations from this scenario occur at the two extreme 
densities. At high density the maximum in the resistivity 
disappears gradually and a purely metallic temperature 
dependence is observed over the entire T-range. At the 
lowest densities one enters the strongly localized regime 
where the presented model can no longer be applied. 
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ii. the effect of energy averaging comes into play and 



FIG. 5. Comparison between experiment and calculations. 
In this calculations we include interface charge and phonon 
scattering. We consider the collisional broadening effect on 
screening with T D = 0.8, 0.6, 0.52, 0.42, 0.37, 0.34T F (from 
top to bottom). In inset the effect of level broadening is 
shown for the density p = 1.95xlO n crn~ 2 . Here only in- 
terface charge scattering is considered. 

We are now able to compare our calculations with the 
experiment. This is illustrated in Fig. a In order to 
achieve quantitative agreement between experiment and 
calculations we additionally included level-broadening ef- 
fects on the polarizability function characterized by the 
parameter T — h/r q = k^To (Tb Dingle temperature, 
marked with filled squares in Fig. nl) Mj. This acts sim- 
ilar to thermal effects in rounding off the sharp corner of 
the polarizability function visible in Fig. pi The influ- 



ence on the resistivity can be seen in the inset in Fig. g: 
it reduces the effect of temperature dependent screening 
and the overall resistance rises. 

The agreement of our simulations with the experiment 
is, even quantitatively, quite good. Only at the lowest 
densities the model fails in describing the experimental 
data, due to the reasons mentioned above. 

In summary we have shown that the non-monotonicity 
in the temperature dependence of the resistivity of p- 
type SiGe samples showing a MIT can be described by 
the temperature dependence of the scattering time td- 
The metallic part stems from the energy averaging of the 
polarizability function II(q,T), the insulating part from 
the energy averaging when performing the integral for 
t(E, T). This interplay is not significantly altered by the 
T-dependence of the chemical potential fi(T). Finally we 
compare our calculations with experimental curves and 
find quantitatively good agreement. 
Financial support from ETH Zurich and the Schweiz- 
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